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Abstract. 

The exceptionally strong coupling realizable between superconducting qubits and 
photons stored in an on-chip microwave resonator allows for the detailed study of 
matter-light interactions in the realm of circuit quantum electrodynamics (QED). 
Here we investigate the resonant interaction between a single transmon-type multilevel 
artificial atom and weak thermal and coherent fields. We explore up to three photon 
dressed states of the coupled system in a linear response heterodyne transmission 
measurement. The results are in good quantitative agreement with a generalized 
Jaynes-Cummings model. Our data indicates that the role of thermal fields in resonant 
cavity QED can be studied in detail using superconducting circuits. 

PACS numbers: 42.50.Ct,42.50.Pq,03.67.Lx,85.35.Gv 
1. Introduction 

In cavity quantum electrodynamics jH [21 [31 H] the fundamental interaction of matter 
and light is studied in a well controlled environment. On the level of individual 
quanta this interaction is governed by the Jaynes-Cummings Hamiltonian [5]. In 
early experiments the quantization of the electromagnetic field was observed in cavity 
QED with Rydberg atoms by measurements of collapse and revival [6]. Similarly, 
the \fn scaling of the atom/photon coupling strength with the number of photons n 
has been observed in the time domain by measuring n photon Rabi oscillations using 
coherent states [7j and Fock states [El |9]. In the frequency domain the y/n scaling can 
be extracted from spectroscopic cavity transmission measurements. Initial attempts 
employing pump and probe spectroscopy with alkali atoms pU] were inconclusive. In a 
recent experiment however the two-photon vacuum Rabi resonance was resolved using 
high power nonlinear spectroscopy [TT]. At the same time, we observed the quantum 
nonlinearity by measuring the spectrum of two photons and one artificial atom in circuit 
QED [12] . In these measurements the originally proposed pump and probe spectroscopy 
scheme [10J was used. Similarly the n photon Rabi mode splitting was studied using 
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multi-photon transitions up to n = 5 [13J. In the dispersive regime the photon number 
splitting of spectroscopic lines [H] provides similar evidence for the quantization of 
microwave radiation in circuit QED. Further experimental progress in the time domain 
has enabled the preparation and detection of coherent states [15] , photon number states 
[IB"] and arbitrary superpositions of photon number states [T7] in circuits. Using Rydberg 
atoms quantum jumps of light have been observed [18], [19] and Wigner functions of Fock 
states have been reconstructed [20J. These experiments demonstrate the quantum nature 
of light by measuring the nonlinear yfn scaling of the dipole coupling strength with the 
discrete number of photons n in the Jaynes-Cummings model [2TJ . 

In this letter we extend our earlier spectroscopic study of the y/n scaling [T2] in 
the circuit QED architecture [22] 123]. We present a measurement of the photon/atom 
energy spectrum up to n = 3 excitations using both thermal and coherent fields. A solid 
state qubit with a large effective dipole moment plays the role of the atom. It is coupled 
to weak microwave radiation fields in a coplanar waveguide resonator which provides a 
very large electric field strength per photon. The coupling strength of a few hundred 
MHz [21] , much larger than in implementations with natural atoms, can freely be chosen 
over a wide range of values and remains fixed in time. The qubit level spectrum is in 
situ tunable by applying flux to the circuit and the internal degrees of freedom can be 
prepared and manipulated with microwave fields. 

Similar solid state circuit implementations of cavity QED have enabled a remarkable 
number of novel quantum optics [23] |25l M, d 13 M, W\ US [131 H3 [291 [30] and quantum 
computation [3H [32J, [33] experiments. 

2. Sample and experimental setup 

For the experiments presented here we use a transmon type qubit [311 135] . which 
is a charge-insensitive superconducting qubit design derived from the Cooper pair 
box [36], as the artificial atom. Its transition frequency is given by Ug^/2-K ~ 
■sJ%E c Ej(§)-E c with the single electron charging energy Ec/h ~ 0.232 GHz, the flux 
controlled Josephson energy Ej{<&) = E Jmax | cos (7r<^> / ^o) I an d E J)max /h pa 35.1GHz, 
as determined by spectroscopic measurements. The two characteristic energies Ej 
and Eq define the full level spectrum of the qubit where the eigenenergy of level I 
is approximately given as E x ~ -Ej + y/8E c Ej{l + 1/2) - E C /12{61 2 + 61 + 3) [M]. The 
cavity is realized as a coplanar resonator with bare resonance frequency v r pa 6.44 GHz 
and photon decay rate k/2h pa 1.6 MHz. Details of the resonator design and fabrication 
can be found in Ref . [37] . Optical microscope images of the sample are shown in Fig. [T^,. 
A simplified electrical circuit diagram of the setup is shown in Fig. [TJo. 

3. Generalized Jaynes-Cummings model 

The transmon qubit is a superconducting circuit with a nonlinear energy level spectrum. 
In many experiments the nonlinearity, which can be adjusted by circuit design and 
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Figure 1. Sample and experimental setup, (a) Optical images of the superconducting 
coplanar waveguide resonator (top) with the transmon type superconducting qubit 
embedded at the position indicated. The qubit of dimensions 300 x 30 /im 2 is shown 
on the bottom, (b) Simplified circuit diagram of the experimental setup. The qubit is 
capacitively coupled to the resonator through C g and flux tuned with a current biased 
(/) external miniature superconducting coil. The resonator, represented by a parallel 
LC circuit, is coupled to input and output transmission lines via the capacitors C m and 
C*out . A quasi thermal field can be generated by modulating a microwave carrier tone 
v r with a large bandwidth pseudo-random noise spectrum S^hfte ■ The noise signal is 
filtered, attenuated and added to pump (v g o,i±, vi±,2±) and probe (V r f) tones. The 
resulting microwave fields are further attenuated and applied to the resonator input. 
After amplification with an ultra low noise amplifier the transmitted probe signal v T f 
is downconverted with a local oscillator (LO) and digitized with an analog to digital 
converter (ADC). 

fabrication, is sufficient to correctly model it as a two level system. For the experimental 
results presented in this work it is however essential to treat the qubit as a multilevel 
system taking into account the coupling of all relevant transmon levels to the cavity 
photons. The physics of this multilevel artificial atom strongly coupled to a single mode 
of the electromagnetic field is described by a generalized Jaynes-Cummings model with 
the Hamiltonian 

H = f\w r a ] a+ ^2 (fi^i °hi (p\-i,i & + o^-i.Oj ■ (1) 

l=e,f,h,... 

Here, hwi is the energy of the Z'th excited state \l) of the multilevel artificial atom, u T is 
the frequency of the resonator field and gi-xj, is the coupling strength of the transition 
I — 1 — > I and one photon, a) and a are the raising and lowering operators acting on 
the field with photon number n and &ij = are the corresponding operators acting 
on the qubit states. In Fig. [2j a sketch of the energy level diagram of the resonantly 
coupled qubit-resonator system [y r = v gfi ) is shown for up to three photons n = 0, 1, 2, 3 
and the first four transmon levels / = g,e, f and h. 

Considering the first two levels of the artificial atom at zero detuning (A = 
v r — ^ g ,e = 0) the eigenstates of the coupled system are the symmetric (\g, n) + 
\e,n — l))/y/2 = | n+) and antisymmetric (\g,n) — |e, n — l))/y/2 = \n—) qubit-photon 
superposition states, see Fig. [2j For n = 1, the coupled one photon one atom eigenstates 
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Figure 2. Sketch of the energy level diagram of a resonant (v r — v g . e ) cavity QED 
system. The uncoupled qubit states \g), \e), \f) and \h) (from left to right) and the 
photon states |0), |1), |2), |3), ... (from bottom to top) are shown with black and red 
solid lines. The dipole coupled dressed states are shown in blue and the shifts due to 
the \ f, 0) and \h,0) levels are indicated by red dashed lines. 



are split due to the dipole interaction [38j EH]- These states have been observed in 
hallmark experiments demonstrating the strong coupling of a single atom and a single 
photon both spectroscopically as a vacuum Rabi mode splitting with on average a 
single atom [ID] and also without averaging using alkali atoms [U], in superconducting 
circuits [23] and semiconducting quantum dot systems [121 S3], or in time resolved 
measurements as vacuum Rabi oscillations at frequency 2g gfi with Rydberg atoms [7J [8] 
and superconducting circuits [T5J [TO]. At least in principle however the observation 
of the first doublet |1±) alone, could also be explained as normal mode splitting in a 
classical linear dispersion theory [44J. The Jaynes-Cummings model however predicts 
a characteristic nonlinear scaling of this frequency as ^/n2g gfi with the number of 
excitations n in the system. In the general multilevel case, the higher energy atomic 
levels renormalize the dipole coupled dressed state energies. This causes frequency shifts 
in the excitation spectrum as indicated in Fig. [2] and the simple y/n scaling is slightly 
modified. 

4. Vacuum Rabi mode splitting in the presence of a thermal field 

In this section we investigate vacuum Rabi resonances both in the presence of a weak 
thermal background field and also in the presence of externally applied quasi-thermal 
fields. To our knowledge this is the first experiment where thermally excited multi- 
photon dressed states are studied systematically in a cavity QED system. 

Experimentally, the coupled circuit QED system is prepared in its ground state 
\g, 0) by cooling it to temperatures below < 20 mK in a dilution refrigerator. Measuring 
the cavity transmission spectrum T in the anti-crossing region of the qubit transition 
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frequency v gfi and the cavity frequency v r yields transmission maxima at frequencies 
corresponding to transitions from \g, 0) to the first doublet |1±) of the Jaynes-Cummings 
ladder, see Fig. |3^. On resonance, see Fig. [3]d, we extract the coupling strength of 
g 5) e/27r = 133 MHz. This is a clear indication that the strong coupling limit g 9je 3> k, 7, 
with photon decay and qubit decoherence rates k/27t,7/27t ~ 1 MHz, is realized. Solid 
lines in Fig. [3^, (and Fig. are numerically calculated dressed state frequencies, see 
Fig. ^jp, with the qubit and resonator parameters as stated above. For the calculation, 
the qubit Hamiltonian is diagonalized exactly in the charge basis. The qubit states 
\g) and |e) and the flux dependent coupling constant g g ^ e are then incorporated in the 
Jaynes-Cummings Hamiltonian Eq. ([!]). Its numerical diagonalization yields the dressed 
states of the coupled system without any fit parameters. 

In addition to the expected spectral lines corresponding to the transition from the 
ground state \g,0) to the first doublet states |1±), we observe three lines with very 
low intensities, see Figs. [3^i and b. These additional transitions are visible because 
the system is excited by a small thermal background field with a cavity photon number 
distribution given by the Bose-Einstein distribution. This thermal field is a consequence 
of incomplete thermalization of the room temperature black-body radiation at the input 




Figure 3. Vacuum Rabi mode splitting with a weak coherent probe tone, 
(a) Measured resonator transmission spectra versus external flux $ close to degeneracy. 
Black indicates low and white high transmission T. The solid lines show dressed 
state energies as obtained numerically and the dashed lines indicate the bare resonator 
frequency v r and the qubit transition frequency Vg } e- (b) Resonator transmission T 
at degeneracy as indicated with arrows in panel (a) , with a line fit to four Lorentzians 
in red. (c) Corresponding energy level spectrum (similar to Fig. |2| with observed 
transitions indicated by arrows and small thermal population indicated with gray 
circles. 
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and output ports of the resonator. A quantitative analysis taking into account the two 
photon states |2±) and the presence of the higher energy qubit levels / and h in Eq. (flj) 
yields the transition frequencies indicated by yellow and red solid lines in Fig. [3j For 
this analysis the coupling constants g e j = 184 MHz and g^h = 221 MHz of higher 
energy qubit levels to the cavity mode, obtained from exact diagonalization of the qubit 
Hamiltonian, have been included. We thus identify two of the additional spectral lines 
as transitions between the first |1±) and second |2±) doublet states of the resonant 
Jaynes Cummings ladder, see Fig. ^jp. The lowest frequency additional spectral line 
corresponds to a transition from the antisymmetric doublet state with one photon |1— ) 
to the qubit / level without a photon \f, 0). The details of the thermally excited 
transmission spectrum can be used as a sensitive probe for the cavity field temperature. 
Analyzing the amplitudes of the Rabi splitting spectrum with a quantitative master 
equation model p3], leads to an estimated photon temperature of T r ~ 0.2 K which 
corresponds to a relatively high mean thermal occupation number of n t h — 0.3 photons 
for the data presented in Fig. 3a. Careful filtering and thermalization at the input and 
output ports results in a typical cavity photon temperature of < 90 mK and < 54 mK 
(n t h < 0.03 and n th < 0.003) as reported in Refs. [2Z] and [T5] . 

In order to access the three photon doublet states |3±) of the coupled multi-photon 
multilevel-atom system we use externally applied broadband quasi-thermal fields. In this 
new approach the dressed eigenstates are populated according to a thermal distribution 
depending on the chosen thermal field temperature. This allows to investigate the 
flux dependence of all resolvable spectral lines for a given effective resonator mode 
temperature T T in a single experimental run. The spectrum of a one dimensional black 
body such as the considered cavity is given as Sid(u) = /iz//[exp (/«///c B T) — 1]. At a 
temperature of T r > hu r /kB — 300 mK this energy spectrum is flat with a deviation of 
< 5% within a 500 MHz band centered at the cavity frequency v r ~ 6.44 GHz. It is 
therefore a very good approximation to make use of a white noise spectrum in the narrow 
frequency band of the experimentally investigated transition frequencies centered at u r , 
in order to generate a quasi-thermal field of temperatures T r > 300 mK and populate 
the considered cavity mode with thermal photons. 

In order to generate such a spectrum a carrier microwave tone at frequency v r is 
modulated with a low frequency large bandwidth quasi-random noise spectrum S^hj^g 
using a mixer, see Fig. [T]d. This approximately yields a microwave frequency white 
noise spectrum with a bandwidth of 500 MHz centered symmetrically around the cavity 
frequency v T . Using tunable attenuators, we can adjust the noise power spectral density 
over a wide range of values. For this experiment we adjust it such that the thermal 
population of the cavity mode is on the order of n t h ~ 0.9 corresponding to a temperature 
of T r ~ 0.4 K. This noise spectrum constitutes a sufficient approximation of a black body 
thermal noise source for the considered ID cavity mode, temperature and frequency. At 
the same time, the chosen mean thermal population n t h ~ 0.9 allows to observe all 
allowed transitions between the ground state \g) and the three photon doublet states 
|3±). 
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In the presence of the thermal field, we probe the cavity transmission spectrum 
as a function of flux in the anticrossing region, see Fig. |4^, with a weak coherent 
probe tone z/ r f in the linear response limit. In this limit the weak probe tone is only 
a small perturbation to the field and no multi-photon transitions are induced. In this 
measurement we resolve all allowed transitions between the thermally occupied dipole 
coupled states in the generalized Jaynes-Cummings ladder. The solid lines are again the 
calculated dressed state transition energies which agree well with the observed spectral 
lines. In Fig. |2}d, a cavity transmission measurement at flux $/$o = 0.25 is shown. We 
identify 8 allowed transitions, compare with Fig. [4)3. It follows that the states |1±), 
|2±) and also \ f, 0) are thermally populated. 

In the two-level-atom approximation, transitions between symmetric and 
antisymmetric doublet states are forbidden at degeneracy. In the generalized Jaynes- 
Cummings model the dressed state transition matrix elements are renormalized due to 
higher qubit levels. Numerical diagonalization however shows that the matrix elements 
squared, which are related to the amplitude of the expected spectral lines, are 140 
(6) times smaller for the symmetry changing transitions |1 — ) — > |2+) (|1+) — > |2— )) 
than for the observed symmetry preserving transitions |1 — ) — > |2— ) (|1+) — > |2+)) 
at degeneracy. Similarly, for the transitions |2— ) — > |3+) (|2+) — > |3— )) the 
matrix elements squared are 235 (16) times smaller than the measured transitions 
|2— ) — > |3— ) (|2+) — > |3+)) at degeneracy. Therefore transitions between symmetric 
and antisymmetric doublet states are not resolved in our experiment. Symmetry 
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Figure 4. Vacuum Rabi mode splitting in the presence of a thermal field, (a) Cavity 
transmission T as in Fig. [3] with an additional quasi-thermal field of temperature 
T r ~ 0.4 K applied to the resonator input populating the |1±), |2±) and |/, 0) 
states, (b) Transmission spectrum at Q/$o = 0.25, indicated by arrows in (a), 
(c) Corresponding energy level spectrum (similar to Fig. [2]) with observed transitions 
indicated by arrows and induced thermal population indicated with gray circles. 
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changing transitions populating the antisymmetric states \2— ) and |3— ) have larger 
matrix elements than symmetry changing transitions populating the symmetric states 
|2+) and |3+) because the former are closer in frequency to the qubit levels \ f, 0) and 
\h, 0). Similarly, the matrix element for the transition |1— ) — > \ f,0) is 34 times larger 
than for the transition |1+) — > \f,0) at degeneracy. The latter is therefore also not 
observed in the experimental data. In addition to the transition |1 — ) — > |/, 0), also 
seen in the data presented in Fig. [3j we observe a transmission line which corresponds 
to the transition \ f, 0) — > \f, 1), see Fig. |4Jd. A numerical calculation shows that the 
matrix element is 5 times larger at degeneracy than |2— ) — > \f,l) and 7 times larger 
than \f, 0) — > \h,0) which in principle could also have been observed. All transitions 
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Figure 5. Measured energy level diagram of the three-photon / artificial atom system. 
Measured dressed state energies (blue dots) reconstructed from extracted transition 
frequencies from data in Fig. [4] compared to calculated uncoupled cavity and qubit 
levels (dashed black lines), the calculated dressed state energies in the qubit two- level 
approximation (dotted red lines) and to the corresponding calculation including four 
qubit levels (solid red lines). 
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observed in the experimental data are in qualitative agreement with the calculated 
matrix elements stated above. 

In Fig. [5] the complete measured level spectrum of the bound photon/atom system 
up to the third excitation is shown. To calculate the absolute energies of the levels (blue 
dots) we extract the transition frequencies from data presented in Fig.|4]with Lorentzian 
line fits and add them accordingly. For the first doublet states |1±) we find excellent 
agreement with both a simple two-level atom Jaynes-Cummings model (dotted red lines) 
as well as the generalized multilevel Jaynes-Cummings model (solid red lines). In the 
case of the second |2±) and third doublet states |3±) we find considerable frequency 
shifts with regard to the two level model (compare dotted and solid red lines) but 
excellent agreement with the generalized model taking account the additional qubit 
levels. Furthermore it can be seen in Fig. [5] that the negative anharmonicity of the 
transmon qubit, together with the strong dipole coupling, causes large frequency shifts 
of the antisymmetric dressed levels |2— ) and |3— ) since they are closer in frequency 
to the qubit levels \f, 0) and |/, 1), \h, 0). This leads to a small reduction of the ^fn 
nonlinearity which is in agreement with the numerical results. 

5. Vacuum Rabi mode splitting with two pump and one probe tone 

In order to probe the excitation spectrum of the two and three photon doublet states we 
can also follow the pump and probe spectroscopy scheme similar to the one presented in 
Refs. [TOl [12], where one of the first doublet states |1±) is coherently pumped and the 
transition to the states |2±) is probed. This technique avoids large intra-cavity photon 
numbers which are needed in high-drive and elevated temperature experiments. In 
analogy to the previous section we wait for the system to equilibrate with the cold 
environment to prepare the ground state |g,0). The qubit is then tuned close to 
degeneracy where $/$ ~ 0.25. We weakly probe the resonator transmission spectrum 
as shown in Fig. [6] (blue lines). In a second step we apply a pump tone at frequency 
v gO,i- (^0,1+) occupying the dressed state |1— ) (|1+)) and probe the system again, 
see Fig. [6^l (b) yellow line. Clearly the transitions |1±) — > |2±) become visible at the 
calculated eigenenergy which is indicated with yellow vertical arrows in Figs. [6^,, b and 
c. We also note that the state \ f, 0) is populated by the probe tone via the transition 
|1— ) — > \ f, 0), see red arrows in Figs. [6^1 and c. In a last step we apply two pump tones 
at frequencies z^,o,i- an d ^1-,2-j see Fig- [6] a (green line), or at frequencies v g o,i+ and 
z/i +;2+ , see Fig. [6]b (green line) respectively. The three-photon one-qubit dressed state 
transitions |2±) — > |3±) become visible in the spectrum, see green vertical arrows. At 
the same time, transitions from the ground state are found to saturate considerably 
when the pump tones are turned on, compare the amplitudes of the spectral lines at 
the frequency indicated by the left blue arrow in Fig. [6] a, or similarly by the right blue 
arrow in figure Fig. |6]b. This is expected since the occupation probability of the ground 
state is reduced and the transition starts to become saturated when the pump tones are 
turned on. 
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Figure 6. Vacuum Rabi mode splitting with one probe tone and zero, one or two 
coherent pump tones at flux Q/Qq ~ 0.25 close to degeneracy, (a) Measured resonator 
transmission spectra T/T max without a pump tone (blue) with one pump tone t'go.i- 
(ycllow) and with two pump tones fgo.i- & n d V\-^.- (green) applied. The spectra are 
offset by 0, 0.02 and 0.03 T/T max and the boxed area is scaled in amplitude by a factor 
of 3 for better visibility. Vertical arrows indicate numerically calculated transition 
frequencies, (b) Similar measurement of resonator transmission T for the case of 
no (blue), |1+) (yellow) and both |1+) and |2+) symmetric dressed states pumped 
coherently, (c) Energy level spectrum (similar to Fig. [2| with relevant transitions 
indicated by arrows. 



Again, the observed transition frequencies are in good agreement with the calculated 
dressed state transition energies indicated by vertical arrows in Figs. [6^, and b. Similar to 
the experiment described in the last section, additional spectral lines with low intensity, 
see Fig. |6^i and b blue lines, occur because of a small probability of occupation of the 
first doublet due to the residual thermal field. In comparison to the data shown in Fig. [4] 
no transition \f, 0) — > \ f, 1) is observed because the level \ f, 0) is neither thermally, nor 
coherently populated here. 

6. Summary 

We extended our previous work [12J by introducing thermal fields to populate the 
dressed eigenstates in a resonant cavity QED system. In addition to the one and two 
photon/atom superposition states we report a measurement of the three photon doublet 
using both thermal and coherent fields. The results are in good agreement with a 
generalized multilevel-atom Jaynes-Cummings Hamiltonian without any fit parameters. 
Related results have been reported in an atomic system [TT] and in circuit QED |13j . 
Similar effects may be interesting to approach also in semiconducting cavity QED 
systems. It has been shown that cavity QED with superconducting circuits can be a 
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sensitive probe for thermal fields. A more detailed quantitative analysis of the thermally 
excited vacuum Rabi spectra could be of interest in the context of environmentally 
induced dissipation and decoherence, thermal field sensing and the cross-over from the 
quantum to the classical regime of cavity QED. 
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